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EXPERIMENTALLY DETERMINED PENDULUM ANALOGY OF LIQUID 
SLOSHING IN SPHERICAL AND OBLATE-SPHEROIDAL TANKS 


by Irving E. Sumner 
Lewis Research Center 

SUMMARY 

An experimental investigation was conducted to determine the general liquid- sloshing 
characteristics (fundamental frequencies, horizontal or side slosh forces, and damping 
ratios) as well as quantities for a pendulum analogy that would effectively represent the 
fundamental mode of liquid sloshing in unbaffled oblate- spheroidal and spherical tanks 
over a range of liquid depths. Tanks having a diameter of 32. 0 inches were used. Ver- 
tical and horizontal slosh forces were measured to determine several of the pendulum- 
analogy parameters. These parameters included (1) the pendulum mass, (2) the length 
of the pendulum arm, (3) the hinge-point location of the pendulum arm, (4) the maximum 
angles through which the pendulum can oscillate, and (5) the fixed mass. 

The experimental results are presented in terms of dimensionless parameters that 
are, in general, independent of tank size, imposed longitudinal acceleration, and density 
and viscosity of the contained liquid. The experimentally determined fundamental- 
frequency and pendulum- analogy parameters are compared, whenever possible, with 
previously obtained analytical results. 


INTRODUCTION 

Propellant sloshing is a potential source of disturbance critical to the stability and/or 
the structural integrity of space vehicles containing relatively large masses of liquid 
propellants, particularly if the propellant oscillations are coupled with either the attitude- 
stabilization-control frequency or the fundamental bodybending frequency of the vehicle. 
Considerable difficulty is found, however, when an attempt is made to include the hydro- 
dynamic equations of liquid motion in a computer simulation of the vehicle motion for a 
stability analysis. Often, the solutions of the hydrodynamic equations do not exist in a 
form that is readily adapted to computer computation. Therefore, it is desirable to rep- 
resent propellant sloshing as a mechanical analogy (e. g. , either a pendulum or a spring- 
mass system) for which the linear, second-order differential equations of motion are 



readily adapted to digital or analog computer simulations. 

A pendulum model used to simulate effectively the fundamental mode of liquid slosh- 
ing generally consists of a fictitious pendulum mass suspended from a given point in each 
propellant tank to represent the sloshing liquid mass and a fictitious fixed mass in each 
tank to represent the nonsloshing mass. Generally, it is necessary to simulate only the 
fundamental mode of liquid sloshing in the mechanical analogy since the higher mode 
natural frequencies are normally much greater than the attitude-control frequencies and 
since the side forces produced at the higher modes are small. The quantities needed to 
describe the pendulum model of liquid sloshing are (1) the pendulum mass, (2) the length 
of the pendulum arm, (3) the hinge-point location of the pendulum arm, (4) the maximum 
angles through which the pendulum can oscillate, (5) the fixed mass, and (6) the centroid 
location of the fixed mass. Analytically determined values of quantities (1), (2), (3), 
and (6) are presented in reference 1 for oblate-spheroidal tanks. There have been no 
known investigations initiated, however, to determine experimentally the values of these 
quantities. Therefore, an experimental investigation was conducted at the NASA Lewis 
Research Center to determine the general liquid- sloshing characteristics (fundamental 
frequencies, horizontal or side slosh forces, and damping ratios) as well as the six pre- 
viously mentioned quantities necessary to represent effectively the fundamental mode of 
liquid sloshing as a pendulum analogy. 

The two tank configurations investigated were (1) an oblate- spheroidal tank having 
a circular diameter of 32. 0 inches and a minor axis of 22. 86 inches and (2) a spherical 
tank having a diameter of 32. 0 inches. Experimental data were obtained over a range of 
liquid-depth ratios and excitation amplitudes for each tank configuration. The contained 
liquid was water. The experimental results are presented in terms of dimensionless 
parameters. The pendulum- analogy parameters are compared, whenever possible, with 
analytical results presented in reference 1. 

APPARATUS AND INSTRUMENTATION 

An oblate- spheroidal tank having a circular diameter of 32.0 inches and a minor axis 
of 22. 86 inches and a spherical tank having a diameter of 32. 0 inches were fabricated from 
clear plastic. The contained sloshing liquid was water in all cases. 

The experimental test facility, which is nearly identical to that described in refer- 
ences 2 and 3, is shown in figure 1. Each of the two tank configurations was mounted on a 
test bed that was suspended from a frame through three vertically oriented load cells and 
one horizontally oriented load cell. The frame was suspended from overhead crossbeams 
and was free to oscillate in one direction in the horizontal plane. The driving force was 
provided by a hydraulic piston and cylinder actuated by an electrically controlled 
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A, Preloading bracket 

B, 10,000-Pound preloaded horizontal load cell 

C, 500-Pound vertical load cel! 

D, Support members from overhead crossbeams to frame 

E, 32 Inch-diameter spherical tank 

F, Stationary support 

G, Hydraulic cylinder 

H, Test bed 

I, Frame 

1, 2, 3 Vertical toad cells 



(a) Schematic view. 



Figure 1. - Experimental slosh-force test facility. 
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servovalve. The excitation amplitude could be varied from 0 to 1 inch, and the excita- 
tion frequency could be varied from 0 to 20 cps. A sinusoidal excitation waveform was 
used for this investigation. The electric and hydraulic control circuits for the driving 
mechanism were designed to enable the oscillatory motion of the frame, test bed, and 
t ^nk to be ’’quick-stopped” at a point of zero velocity during any given cycle of oscilla- 
tion so that only the residual forces resulting from the liquid sloshing could be measured. 

The horizontal load cell was a piezoelectric quartz crystal that had been preloaded 
in compression. The vertical load cells were semiconductor-type transducers. The 
residual horizontal and vertical forces resulting from the liquid sloshing were sensed by 
the horizontal and vertical load cells, respectively, and the signals were displayed as 
continuously recorded oscillograph traces. The maximum error in obtaining the slosh 
forces from the oscillograph traces was approximately ±2. 5 percent of the actual value. 


PROCEDURE AND DATA REDUCTION 
General Liquid-Sloshing Characteristics 

Each tank was oscillated sinusoidally at a preselected excitation frequency (0. 66 to 
1. 93 cps) that encompassed the range of fundamental frequencies of liquid motion over 
the range of liquid depths investigated and at an arbitrarily selected amplitude (0. 010 to 
0. 200 in. ). At each liquid depth investigated, the excitation frequency chosen was equal 
to the fundamental frequency of oscillation of the contained liquid since this provided the 
maximum wave height of the liquid surface and, therefore, the maximum slosh force at 
a given excitation amplitude (refs. 3 and 4). These maximum slosh forces are herein- 
after referred to as the first-mode slosh forces. The maximum wave height of the liq- 
uid surface was established (1) for a steady- state oscillation of the liquid surface at low 
liquid depths and (2) for the point where the relatively flat liquid surface began to break 
up and shower liquid throughout the tank and before any rotary or swirling motion of the 
liquid could take place at the higher liquid depths. When the fundamental or first natural- 
mode waveform had built up to its maximum height on the tank wall, the oscillatory mo- 
tion of the tank was quick- stopped, and the residual slosh forces were recorded. 

The fundamental frequencies of the liquid oscillations were determined for the first 
several slosh-force cycles occurring immediately after the qui ck s top and are presented 
in the form of the fundamental-frequency parameter r] e = Vr/g. (Symbols are de- 
fined in appendix A. ) The values of the first- mode horizontal slosh forces were deter- 
mined from the first force peak occurring immediately after the quick stop; the hori- 
zontal forces are presented in the form of the first- mode slosh force parameter 

A = F /pgaD . The first- mode damping ratios were calculated by averaging the first 
s 
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two successive values of the logarithmic decrement 5 = In [(F s ) n /(F s ) n+ jJ on each oscil- 
ograph trace, where (F ) was the peak force on one slosh cycle and (F ) 1 was the 

peak force on the succeeding cycle. 


Quantities for Pendulum Analogy 

Pendulum mass . - In order to determine experimentally the pendulum mass (or ef- 
fective sloshing mass), it was necessary to oscillate each tank at an excitation frequency 
much less than the fundamental frequency of the contained liquid such that, for steady- 
state conditions, the liquid would oscillate at exactly the excitation frequency. Excita- 
tion frequencies were varied from 40 to 70 percent of the fundamental frequencies over 
the range of liquid- depth ratios investigated, and excitation amplitudes were varied from 
0. 200 to 0. 900 inch. Specific values of the excitation frequency and amplitude at each 
liquid-depth ratio were chosen such that the wave height of the liquid surface was as 
large as possible while the liquid surface remained relatively flat. Thus, the accuracy 
of measuring the relatively large slosh forces present was increased. Once the wave 
height reached a steady-state value, the oscillatory motion of the tank was quick- stopped, 
and the horizontal slosh forces were determined from the first force peak occurring 
immediately thereafter. The method of calculating the pendulum mass is presented in 
appendix B. 

Hinge -point location . - The excitation frequency was equal to the fundamental fre- 
quency of the contained liquid. The excitation amplitude was varied from 0. 005 to 0. 050 
inch depending on the liquid- depth ratio. The wave height was allowed to increase toward 
a maximum value with the following restrictions: (1) the liquid surface would remain 
relatively flat with no splashing of the contained liquid and (2) the liquid would be oscil- 
lating in exactly the same direction as the driving force (no swirl or rotary motion of the 
liquid surface) so that the slosh-force traces of vertical load cells 1 and 2 (fig. 1(a), 
p. 3) would be in phase and equal in magnitude. The oscillatory motion of the tank was 
then quick-stopped, and the residual forces resulting from the liquid sloshing were re- 
corded. 

A typical slosh-force oscillograph trace is shown in figure 2. The oscillograph 
trace labeled Tank position indicates (1) the initial oscillatory motion of the frame, test 
bed, and tank and (2) the point at which the oscillatory motion was quick-stopped. The 
data trace labeled F c indicates the slosh-force measurements obtained from the hori- 
zontal load cell, while the traces labeled F., Fg, and Fg indicate the force measure- 
ments obtained from vertical load cells 1, 2, and 3, respectively. (The position of the 
vertical load cells is noted in fig. 1. ) The horizontal and vertical slosh forces were de- 
termined from the first force peaks occurring immediately after the tank was quick- 
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Figure 2. - Typical slosh-force oscillograph trace for spherical tank. Liquid-depth 
ratio, 0.50; excitation-amplitude parameter, 0.00094. 


stopped, and the moment about the center of the tank was calculated. The method of ob- 
taining the hinge-point location is given in appendix B. 

Liquid-slosh angle. - The angle from the horizontal at which the liquid surface 
oscillated y ^ (fig. 3) was measured visually by a protractor-type device when the liquid 
surface had attained its maximum height on the tank wall for different values of 
excitation- amplitude parameter and liquid- depth ratio. The accuracy of this device was 
better than ±2°. 


EXPERIMENTAL RESULTS 
General Liquid-Sloshing Characteristics 

Fundamental frequency . - By using an "equivalent cylindrical tank" method (ref. 5), 
the fundamental frequencies of oscillation of the liquid contained in the oblate- spheroidal 
and spherical tanks were calculated from the equation 



Figure 3. - Definition of liquid slosh angle. 



( 1 ) 


The experimental values of the correction factor a ob- 
tained in references 2 and 5 to provide closer agreement 
between the calculated and experimental values of the 
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Figure 4. - Ratio of experimental to calculated fundamental frequency for oblate- 
spheroidal tank (refs. 2 and 5). 


fundamental frequency over a wide range of liquid-depth ratios are shown in figure 4. 

The values of the fundamental frequency had also previously been obtained analytically 
from exact solutions of the hydrodynamic equations of liquid motion for a spherical tank 
(ref. 4). The analytical, calculated, and experimental values of the fundamental fre- 
quency are presented in terms of the parameters r ; a = o> a V r /g> V c = <^ c V r /g, and 
rj e = 5 Vr/g, respectively, in figure 5 for the oblate- spheroidal and spherical tanks. 

The fundamental-frequency parameter increased with liquid-depth ratio for the two tank 
configurations investigated. The analytical, calculated, and experimental values were in 
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(a) Oblate-spheriodal tank. 


- r r — | — | — | — j—| — | — | — | close agreement, with less than 4 percent variation 

- - amplitude" j — in all cases. Each data point shown represents an 

1 parameter, average of several values; the maximum deviation 

1 X 0 /D 

„ from the average was ± 1 . 5 percent. 

I o Q. OOoo _i _L 

° .00313 Horizontal slosh forces. - The horizontal slosh 

I — o .00156 — 

a ,00063 forces were obtained for an excitation frequency 

v I 'i 00031 equal to the fundamental frequency of the contained 

Q5 _ liquid at each liquid-depth ratio investigated and are, 

therefore, the maximum slosh forces that could be 

w Lvf \ 'vt — obtained for a given excitation amplitude. These 

[J] / ^ \ ikL slosh forces are presented in figure 6 in terms of the 

03 S, first-mode slosh-force parameter X = Fg/pgaD^ 

„ 6W j /Vvk Hi for the oblate-spheroidal tank and the spherical tank. 

JL .02 Mjy f The maximum slosh forces for a given excitation 

A J / a . I * amplitude were found when each tank was at or very 

■s' ‘ _ ffi J y \ near a half -full condition. The slight dip that ap- 

I nl 111111 IV; peared in the faired curves at liquid- depth ratios of 

f 06i a) Obiate-spheriodai tank. o. 0 q to 0. 75 and for excitation- amplitude parameters 

t — — equal to or less than 0.00156 occurred because the 

■I .05 liquid near the surface tended to separate from the 

| ty\\\ tank wall as the wave height of the liquid surface ap- 

I; .M /t/y® > \ > proched the maximum height that could be obtained. 

i (J \ To maintain a relatively flat liquid surface, it was 

H [ /\ \ ° 'Ml necessary to quick-stop the oscillatory motion of the 

02_ - if] f tank somewhat prematurely, which thereby reduced 

ly A the slosh forces slightly. 

01 — -Ju J — Damping ratios . - The first- mode damping ra- 

- wjys-* lkr tios Gained in the oblate-spheroidal and spherical 

^ — 4 — g *'■' — ^ — j 0 tank configurations are presented in figure 7 for a 
Liquid -depth ratio, h/aD range of excitation amplitudes. At low liquid- depth 

(b) Spherical tank. ratios, the damping ratios increased with a decrease 

i" the excitation amplitude and liquid-depth ratio be- 
cause of the small mass of liquid oscillating over a 
large tank wall area. At high liquid-depth ratios, the damping ratios increased with an 
increase in excitation amplitude and liquid-depth ratio because of (1) the curvature of 
the tank wall and (2) the small portion of liquid that splashed from the tank wall and 
tended to produce ripples on the surface of the liquid. 

Scaling factors . - The fundamental frequencies of liquid oscillations in oblate- 
spheroidal or spherical tanks of any size under any imposed longitudinal acceleration can 
be calculated directly from equation (1) (p. 6), or from figure 5 for geometrically similar 


Liquid-depth ratio, h/aD 
(b) Spherical tank. 

Figure 6. - First-mode slosh -force parameter 
as function of liquid-depth ratio. 
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lb) Spherical tank. 


Figure 7. - First-mode damping ratio as 
function of liquid-depth ratio. 


tank configurations. The fundamental frequencies 
are generally independent of the kinematic viscosity 
of the contained liquid (ref. 3). 

When it is necessary to predict sloshing char- 
acteristics in a full-scale tank from experimental 
tests conducted in a scale- model tank, it is desir- 
able for the values of the viscosity parameter 
v/\ f/gD 3 for the two tanks to be nearly the same. 

But in cases where the full-scale tank is appre- 
ciably larger than the model tank, it is generally 
impossible or impractical to obtain a sloshing liquid 
with a low enough kinematic viscosity for the model 
tank. The slosh-force parameter is virtually inde- 
pendent of tank size and liquid viscosity, however, 
when the tank diameter is large or when the liquid 
viscosity is small (values of the viscosity parameter 
less than 10 ). This is indicated in figure 8 for 

half -full spherical tanks oscillated at an excitation- 
amplitude parameter Xq/D = 0. 01 (ref. 3). The 
value of the viscosity parameter for a full-scale 
tank is generally much less than 10 . If the model- 

tank size and the slosh liquid are chosen such that 
the viscosity parameter of the model is also less 
than 10” 6 , the sloshing characteristics (except for 
the liquid damping) in the full-scale and model tanks 
will be very similar, and the effect of liquid viscosity 
on the slosh-force parameter may be neglected. The 



Figure 8. - First-mode slosh-force parameter as function of viscosity 
parameter for spherical tanks. Liquid-depth ratio, 0.5; excitation- 
amplitude parameter, 0.01. 
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Figure 9. - First-mode damping ratio as function of viscosity parameter 
for spherical tanks. Liquid-depth ratio, 0.5. 



first-mode slosh forces expected to occur in 
geometrically similar tanks can be obtained 
from the values of the force parameter pre- 
sented in figur e 6 (where the viscosity param- 
eter v/ |/Jd3 = 4. 3X10- 7 ). 

It is somewhat more difficult, however, to 
predict from scale- model-tank tests the first- 
mode damping ratios that are expected to occur 
in a full-scale tank when the values of the vis- 
cosity parameter for the two tanks are not 
nearly the same. In this case it is usually nec- 
essary to extrapolate experimental data obtained 
from several model -tank tests where the viscosity parameter has been varied over a sub- 
stantial range. Figure 9 shows the viscous first-mode damping ratios obtained in half- 
full spherical tanks (ref. 3). In general, the damping decreases for any tank configura- 
tion (1) with a decrease in the liquid kinematic viscosity or (2) with an increase in the 
tank size or longitudinal acceleration imposed on the tank. 


Figure 10. - Quantities used to describe pendulum analogy of 
liquid sloshing. 


Pendulum Analogy 

The physical quantities necessary to represent the fundamental -mode liquid-sloshing 
characteristics as a pendulum analogy are shown in figure 10. 

Pendulum (or sloshing) mass . - The effective pendulum mass that would produce 
oscillatory side forces equal to those that resulted from the sloshing liquid was experi- 
mentally determined by using equation (B6). The ratio of the sloshing mass to the total 
liquid mass present at a given liquid-depth ratio nip/mj. is shown in figure 11 for the 
oblate- spheroidal and spherical tank configurations. At low liquid-depth ratios, almost 
all of the liquid mass in each tank acted as the sloshing mass, while at the higher depth 
ratios, only a small portion of the liquid mass acted as the sloshing mass. The 
sloshing- mass ratio was greater for the oblate -spheroidal tank at a given depth ratio. 
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T T 7 1 

Excitation- Excitation - 
frequency amplitude 

ratio, parameter, 



0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

Liquid-depth ratio, h/aD 


Figure 11. - Ratio of pendulum or sloshing mass to total liquid mass in partially filled tank. 



Figure 12. - Ratio of pendulum or sloshing mass to total liquid mass in completely 
filled tank. 

No dependence of the sloshing- mass ratios on either the excitation frequency or the ex- 
citation amplitude was noted. 

The sloshing mass ratio m/m— for a range of liquid depth ratios 0 <C h/ a© < 1. 0 
is presented in figure 12 where the sloshing or pendulum mass m for a given liquid 

Sr 

depth has been arbitrarily divided by a constant, namely the total liquid mass m,p con- 
tained in a completely filled tank. The curves shown here represent values obtained from 
the faired curves presented in figure 11. The maximum sloshing mass was obtained at a 
liquid-depth ratio of 0. 6 for both tank configurations, although the maximum slosh forces 
had been obtained at a depth ratio of 0.5 (fig. 6). 
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Figure 13. - Pendulum-arm-length ratio. 



Pendulum arm length. - The pendulum arm length for each tank configuration was 
calculated directly from the pendulum fundamental-frequency expression Lp = gA> e ; 
the ratio of pendulum arm length to tank height Lp/aD is presented in figure 13. The 
pendulum-arm-length ratio decreased as the depth ratio increased and was smaller for 
the spherical tank at a given liquid-depth ratio. 

Hinge-point location of pendulum arm . - The location of the hinge point of the pendu- 
lum arm was assumed to lie on the vertical axis of symmetry of the tank, and its loca- 
tion on that axis was determined by using equation (B9). The ratio of hinge-point location 
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Maximum liquid-slosh angle, 7 t , or pendulum angle, y p , deg 


to tank height H^/aD is presented in figure 14 for the two tank configurations. The 
hinge-point location moved from near the top to the center of the oblate- spheroidal tank 
as the liquid-depth ratio increased from 0 to 1. 0. The hinge-point location tended to re- 
main at the center of the spherical tank regardless of the liquid-depth ratio. 

Pendulum angle . - The validity of equation (B6) (appendix B), from which the pendu- 
lum mass was determined, depended on the assumption that the pendulum arm remained 
normal to the liquid surface. The angles y^ from the vertical, through which the pen- 
dulum mass must oscillate to produce side forces equal to the first-mode horizontal 
slosh forces, were calculated from equation (B2) for a range of liquid-depth ratios 

(0.2 <h/«D < 0. 8) and excitation am- 
plitudes (0. 010 to 0. 200 in. ) for both 
tank configurations. The slosh angles 
y^ were measured visually for the 
same initial sloshing conditions 
(o>q and Xq) as for the calculated 
values of the pendulum angles y . A 
comparison of the calculated (y ) and 
measured (y^) angles is shown in fig- 
ure 15 for the oblate and spherical 
tanks. Agreement was generally good, 




0 1 Z 3 4 5 6 7x10-3 

Excitation -amplitude parameter, Xq/D 


and the assumption that the pendulum 
arm remained normal to the surface of 
the liquid (appendix B) was thereby 
confirmed. 


Fixed (or nonsloshing) mass. - 


The fixed mass was determined from 


equation (B10) and is presented in the 


form of the fixed- mass ratio m^/m^, 
in figure 16 for the two tank configura- 


tions. The fixed-mass ratio increased 


with an increase in liquid- depth ratio. 
The curves shown here represent val- 
ues obtained from the faired curves 


presented in figure 11. 

Scaling factors . - The length of 
the pendulum arm depends only upon 
the fundamental frequency of the con- 
tained liquid. Therefore, the 
pendulum- arm-length ratio for 


Figure 15. - Comparison of liquid-slosh angle and pendulum angle. 
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geometrically similar tanks is independent of the tank size, the imposed longitudinal 
acceleration, and generally, the kinematic viscosity of the contained liquid. 

The values of the pendulum and fixed- mass ratios, the hinge -point -location ratio, 
and the pendulum angles that are presented should be independent of the tank size and im- 
posed longitudinal acceleration for geometrically similar tanks for values of the vis- 

0 

cosity parameter less than approximately 10 

Comparison of Experimental Pendulum Analogy 
Parameters With Analytical Values 

An analytical solution to the hydrodynamic equations of liquid motion and the equiva- 
lent pendulum-analogy parameters for oblate- spheroidal tanks of varying eccentricity 
(tank height-diameter ratio) were obtained by A. H. Hausrath. The results of this in- 
vestigation are presented in reference 1 for tank height- diameter ratios a of 0. 990, 

0. 752, and 0. 500. A comparison of these analytical results with the experimental re- 
sults obtained in this investigation for tank height-diameter ratios of 1. 00 and 0. 715 
showed the following: 

(1) The analytical and experimental values of the fundamental-frequency parameter 
(fig. 17) agreed fairly well except that (a) some difference in results was noted at liquid- 
depth ratios above 0. 7 and (b) analytical values for a = 0. 752 appeared to be low for 
liquid-depth ratios less than 0. 7. 

(2) Analytical and experimental values of the pendulum-mass ratio m^/m^, (fig. 18) 
tended to reach a maximum at approximately the same liquid- depth ratio. However, the 
experimental values were as much as 56 percent greater at some depth ratios. 



0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

Liquid-depth ratio, h/aD 

Figure 18. - Comparison of experimental pendulum- or s!oshing-mas$ ratio with analytical values. 
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Figure 21. - Analytical values of centroid location of fixed or nonsloshing mass. (Data 
from ref. 1. ) 


(3) Analytical and experimental values of the pendulum-arm-length ratio (fig. 19) 
showed good agreement. 

(4) Analytical and experimental values of the pendulum-arm-hinge-point-location 
ratio (fig. 20) showed good agreement. 

Analytical values only of the centroid location of the fixed mass are presented in fig- 
ure 21 since it was not possible to determine this quantity experimentally. The centroid 
location ratio £ cg /aD generally tended to decrease with an increase in (1) the tank 
height- diameter ratio a and (2) the liquid-depth ratio. The spherical-tank configuration 
would be an exception in that the centroid location of the fixed mass should remain at the 
center of the tank regardless of the liquid-depth ratio. 

SUMMARY OF RESULTS 

An experimental investigation was conducted to determine the general liquid- sloshing 
characteristics and a pendulum analogy that would effectively represent the fundamental 
mode of liquid sloshing in oblate-spheroidal and spherical tanks over a range of liquid- 
depth ratios. The experimental results are presented in terms of dimensionless param- 
eters and are compared, whenever possible, with analytical results. These results may 
be summarized as follows: 

1. The fundamental-frequency parameters increased with an increase in the liquid- 
depth ratio. Values of the analytical, calculated, and experimental fundamental- 
frequency parameters were within 4 percent in all cases. 

2. The maximum first- mode slosh-force parameters obtained for the oblate- 
spheroidal and spherical tanks were found when each tank was at or very near a half-full 
condition. 

3. The first- mode damping ratios for low liquid-depth ratios increased with a de- 
crease in the excitation- amplitude parameter or liquid-depth ratio. At high liquid-depth 
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ratios, the damping ratio increased with an increase in the excitation-amplitude param- 
eter or liquid- depth ratio. 

4. The ratio of pendulum or sloshing mass to total mass in the partially filled tank 
decreased as the liquid-depth ratio increased. The ratio of sloshing mass to total mass 
in the completely filled tank attained maximum values at a liquid-depth ratio of about 0. 6. 
Values of both ratios were slightly greater for the oblate- spheroidal tank at a given 
liquid- depth ratio. 

5. The ratio of pendulum arm length to tank height decreased as the liquid-depth ra- 
tio increased; values of this ratio were greater at a given liquid-depth ratio for the 
oblate- spheroidal tank. 

6. The hinge-point location of the pendulum arm moved from near the top to the cen- 
ter of the oblate-spheroidal tank as the liquid-depth ratio increased, while the hinge- 
point location remained in the center of the spherical tank regardless of the liquid-depth 
ratio. 

7. Close agreement between the calculated values of the pendulum angle and the mea- 
sured slosh angles indicated that the pendulum arm remained normal to the liquid surface 
over a range of liquid- depth ratios for both tank configurations. 

8. The ratio of fixed mass to total mass in the completely filled tank increased as 
the liquid-depth ratio increased. 

9. Experimental values of the fundamental-frequency, pendulum- arm-length, and 
hinge -point- location parameters for both tank configurations showed good agreement with 
analytical values obtained previously. Experimental values of the ratio of sloshing mass 
to total mass in the completely filled tank showed the same characteristics but were 
somewhat greater than the analytical values. Analytical values of the fixed-mass 
centroid-location parameter generally tended to decrease with an increase in (1) the tank 
height- diameter ratio and (2) the liquid-depth ratio. 

10. The values of the fundamental-frequency parameter and the pendulum-arm-length 
ratio are generally applicable to tanks (1) of any size having similar geometric configura- 
tions and (2) under any imposed longitudinal acceleration; the kinematic viscosity of the 
contained liquid generally should not have any effect. All the remaining liquid- sloshing 
and pendulum- analogy parameters presented, with the exception of the damping ratios, 
are believed to be applicable under conditions (1) and (2) as long as the value of the vis- 
cosity parameter is less than approximately 10 Appropriate values of damping in 
large-scale tanks must generally be extrapolated from small-scale tests in geometrically 
similar tanks for a varying viscosity parameter. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 6, 1964. 
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APPENDIX A 
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h/«D 
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eg 


cg,q 


SYMBOLS 


minor axis of oblate - 
spheroidal tank, ft 

circular diameter of tank, ft 

force on tank producing pro- 
pellant sloshing, lb 

horizontal (side) slosh force, 
lb 

vertical acceleration of 
tank, 32. 174 ft/sec ^ 

liquid depth, ft 

liquid depth in flat-bottom 
cylindrical tank of radius r 
containing same liquid vol- 
ume as partially filled 
oblate- spheroidal or spher- 
ical tank, - MhJjfP - h >l, 

3 L orD - h J 

ft 

liquid-depth ratio 

moment of inertia of fixed 
(nonsloshing) mass, slugs/ 
sq ft 

length of pendulum arm, ft 

pendulum-arm-length ratio 

distance from center of tank 
to centroid of fixed mass, ft 

distance from centroid of 
fixed mass to vehicle center 
of rotation 


{ /«D 

eg 7 



H p /«D 
M 


m 


P 


m n 


m. 




q 

r 


t 

x o 

x 0 A> 


centroid location of fixed- 
mass ratio 

distance from center of tank 
to hinge point of pendulum 
arm, ft 

distance from hinge point of 
pendulum arm to vehicle 
center of rotation 

pendulum-arm-hinge-point- 
location ratio 

external moment on tank pro- 
duced by liquid sloshing, 
ft-lb 

pendulum (effective liquid 
sloshing) mass, slugs 

total liquid mass present in 
completely filled tank, slugs 

total liquid mass present in 
partially filled tank, slugs 

fixed (effective liquid non- 
sloshing) mass, slugs 

vehicle center of rotation 

radius of liquid surface of 
partially filled tank, ft 

time, sec 

maximum excitation ampli- 
tude, ft 

excitation- amplitude param- 
eter 
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X 


x 0 

a 


e 


V 


a 






e 


direction of force on tank pro- 
ducing propellant slosh 

excitation amplitude at any 
time t, Xq sin a> O 1 ’ ft 

ratio of tank height to diam- 
eter, b/D for oblate 
spheroid, D/D = 1 for 
sphere 

angle from horizontal through 
which liquid surface oscil- 
lates, radians 

angle from vertical through 
which pendulum oscillates, 
radians 

first- mode damping ratio 
(logarithmic decrement), 

* n L^s^n^^Vn+l] 

first root of JJ (e n ) = 0, 

1. 841 

analytical fundamental- 
frequency parameter, 



fundamental-frequency 
parameter calculated from 
equivalent-cylindrical-tank 

method, u) c j/r/g 

experimental fundamental- 
frequency parameter, 

angular rotation about q, 
radians 


A slosh-force parameter, 

F s /pgaD 3 

v liquid kinematic viscosity, 

sq ft/sec 

v/ yS? viscosity parameter 

£ damping factor, 6 /2 tt 

p liquid mass density, 

slugs/cu ft 

a ratio of experimentally 

determined to calculated 
fundamental frequency 
(refs. 2 and 5) 

oj fundamental frequency of 

liquid oscillation, radians/ 
sec 

u>q excitation frequency of tank, 

radians /sec 

Subscripts: 

a analytically determined 

(ref. 4) 

c calculated by equivalent 

cylindrical-tank method 
(refs. 2 and 5) 

e experimentally determined 

n cycle number, 1, 2, 3, . . . 

Superscripts: 

• first derivative with respect 

to time 

" second derivative with re- 

spect to time 
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APPENDIX B 


PRESENTATION OF LIQUID SLOSHING AS PENDULUM ANALOGY 


Figure 22 shows a pendulum analogy of liquid sloshing for an arbitrary propellant- 
tank configuration symmetrical about its longitudinal axis and fixed in a rigid-body ve- 
hicle or test facility. 

The summation of forces acting on the 
tank in the x- direction that result in liquid 
sloshing may be written as 



T 


P,Q 


‘■eg, q 


F = - m A x A - m A f 9 f m gy (Bl) 
0 0 0 eg, q p &A p ' ’ 


Small angle approximations have been as- 
sumed in all equations. In this investiga- 
tion, the tank and test facility were not sub- 
jected to any rotational, or pitching, motion 
so that 9 = 0 . In addition, the normal test 
procedure was to oscillate sinusoidally the 
tank and the contained liquid at a given ex- 
citation amplitude and frequency until the liquid surface had attained its maximum wave 
height. The motion of the tank was then quick- stopped, after which x Q = 0, and the 
forces resulting only from the liquid sloshing were experimentally measured (F = F ). 
Equation (Bl) may then be written as 


Figure 22. - Pendulum analogy of liquid sloshing. 


F = m gy 
s P P 


(B2) 


Equation (B2) would be useful in calculating the pendulum mass if the pendulum angle 
corresponding to a specific slosh force F g were known. The angle could be de- 
termined directly if the pendulum arm were assumed to remain normal to the liquid sur- 
face (y p = 7f). It was difficult, however, to obtain sufficiently accurate measurements 
of the liquid-slosh angle visually, and a second method of determining the pendulum 
mass was utilized. 

The equation of the pendulum motion may be written as 


m L y + 2m L y + m L co y = -m x. - m (St „ - L )0 
P P P p s e p A p p p e A p p 0 p P,q P 


(B3) 
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If the damping factor 
negligible, and since 


| is small, the damping term ( 2ra p^ e L p r p ) can be considered 
0 = 0 for this investigation, equation (B3) may be rewritten as 




+ u> v ) = -m x_ 
e A p ; p 0 


(B4) 


When the liquid mass was oscillated at an excitation frequency much lower than its fun- 
damental frequency (Wq « o> e ), the liquid oscillated at exactly the excitation frequency 

as long as the oscillatory motion of the tank continued. Since y n was assumed to be 

.... 2 P 

equal to y^, it can also be assumed that y^ = y^ = -o>q y^ sin u>Qt. At the maximum 
wave height of the liquid surface, then, it is assumed that y^ = - w Q>p- By making this 
substitution, equation (B4) becomes 


t / 2 2> 

m L y (-u + u> ) = -m x n 
p p'p' 0 e' p 0 


(B5) 


Equation (B5) must again be rewritten in terms of measurable quantities by making the 
following substitutions: 

2 

(ULp = g/^e (pendulum frequency expression) 

M 2 

(2) Xq = -WqXq (at maximum wave height of liquid surface) 

(3) m = F /g y (substitute into left side only of eq. (B5)) 

p s p 

The resulting equation is then 


m p = 


F 

s 




(B6) 


Equation (B6) is valid for determining the pendulum or effective-sloshing mass as long 
as the damping factor £ is small and, therefore, the damping term in equation (B3) is 
negligible. Consequently, the values of the first force peaks occurring immediately 
after the quick stop are almost the same values as the steady-state slosh forces that 
occurred while the tank was subjected to the oscillatory motion. 

Once the pendulum mass has been determined for a given liquid-depth ratio, the 
pendulum angle y can be calculated for any corresponding value of the horizontal slosh 

r 

force by using equation (B2). 

The hinge-point location of the pendulum arm may be determined from a summation 
of moments about the vehicle center of rotation q: 


0 cg,q 0 


i x n - 


0 + m 0 cg,q' y T “* 


p p,q' 


(B7) 
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The test facility was not subjected to any rotational, or pitching, motion (9 = 0), how- 
ever, and a summation of moments about the center of the tank may be written as 


M = - m oVo +m pV , P 


(B8) 


Since nip = F g /gyp anc * Xq = 0 after the tank has been quick-stopped, the location of the 
hinge point above the center and on the longitudinal axis of the tank may be determined 
from 


f p = M/F g (B9) 

No feasible method was available to determine experimentally the fixed or effective 
nonsloshing mass or its effective centroid location. It may be assumed, however, that 
the fixed mass may be determined from the following equation: 

m 0 = m t ‘ m p (B10) 

The ability of a pendulum analogy to represent the fundamental mode of liquid slosh- 
ing in a cylindrical tank is discussed in reference 6 for two cases: (1) translational 
oscillatory tank motion, and (2) pitching oscillatory tank motion. Additional information 
on the numerical values of the pendulum analogy parameters for spherical and cylindrical 
tanks, as well as the effect of liquid sloshing upon vehicle stability, can be obtained from 
other sources (refs. 7 and 8) as well as from reference 1. It should be noted, however, 
that some differences in numerical values for the fundamental frequency parameter, the 
sloshing or pendulum mass, and the effective centroid location of the fixed mass were 
found to exist between the experimental and/or analytical results presented herein and 
those presented in reference 8. 
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